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PACS. 74.70.Pq - Ruthenates. 
PACS. 74.25.Gz - Optical properties. 

PACS. 74.20.Rp - Pairing symmetries (other than s-wave) 



Abstract. - We compute the optical conductivity of 2D f-wave superconductors and also 
of the multigap model proposed recently by Zhitomirsky and Rice at T = OK in the Born 
limit. The presence of interband impurity scattering was found to play an important role: the 
contributions from the two bands mix up, and new structures are seen in the tunneling density 
of states and in the optical spectrum as well, corresponding to interband transitions. This will 
provide a sensitive test in selecting the competing models for the triplet superconductivity in 
Sr 2 Ru0 4 . 



Introduction. - The discovery of superconductivity in Si'2Ru04 in 1994 has generated 
much interest since it is a perovskite with the same crystal structure as La2Cu04 £Q. From 
analogy with superfluid 3 He Rice and Sigrist |2H3] proposed the spin triplet p-wave supercon- 
ductor with order parameter 

A(k) = dA p (T)(k x ±ik y ) (1) 

Here d is the spin vector parallel to z. 

The triplet pairing has been confirmed by the presence of spontaneous magnetization seen 
by /iSR 0] , which revealed that the superconducting state breaks the time- reversal symmetry, 
and a flat 17 Knight shift seen by NMR 5 showed no change in the spin susceptibility when 
passing through the superconducting transition temperature. However the large density of 
states (N(0)) in the gapless region was somewhat surprising [B], since eq. ^ has the full 
energy gap [7|. 

Furthermore, as sample quality improved, the characteristics of nodal superconductors 
became clearly visible in the specific heat (oc T 2 ) 0, the magnetic penetration depth (oc 
T 2 ) 9], NMR relaxation rate (oc T 3 ) [TU], the thermal conductivity (oc T 2 ) [TT1IT2I and the 
ultrasonic attenuation (TjJJ. Therefore a variety of 2D f-wave models have been proposed 
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|14U15llT5| . Since all these f-wave models have the same quasiparticle density of states, the 
specific heat and the magnetic penetration depth etc. are described equally well by any of 
these 2D f-wave models |16| . On the other hand, the angular dependent magnetothermal 
conductivity data by Izawa et al. |17j can exclude all 2D f-wave models with nodes lying in 
the a- b plane |T5IIT5| . 

Therefore it singles out the most consistent f-wave as 

A(k) = dA f (T)(k x ± ik y ) cos{ck z ). (2) 

The nodal lines are horizontal. Also the pair correlation takes the maximum value for the 
relative separation of the pair at ±c. However in quasi-two dimensional systems with para- 
magnon (i.e. the ferromagnetic fluctuation), the ground state of p-wave superconductors is 
given by eq. Q 3j. This is why Zhitomirsky and Rice (ZR) ^] proposed the multigap model: 
there is p-wave superconductor attached the active 7 band, which induces another f-wave like 
order parameter in the passive a + f3 bands through Cooper pair scattering given by 

A 2 (k) = dA 2 (T)(fc x ± ik y ) cos(ck z /2). (3) 

Since A 2 (k) as given in eq. Ogives the same quasiparticle density of states as a d-wave super- 
conductor |2l)j. it is not difficult to reproduce the specific heat and the magnetic penetration 
depth within the multigap model as was shown recently by Kusunose and Sigrist |2~T]. 

However a simple analysis suggests that A 2 (k) is incompatible with the angular dependent 
magnetothermal conductivity data. In fig. (QJ we show A(k) for p-wave, f-wave and A 2 (k) 
in quasi 2D systems. 




p-wave |A2(k)| f-wave 



Fig. 1 - The absolute value of the p-wave (left), A2(k) (middle) and f-wave gap is shown. Note the 
presence of horizontal nodes in the latter two. 

We note here that v c = (i.e. the Fermi velocity in the z direction) on the nodal lines of 
A 2 (k) which gives rise to a large cos(20) term in the angular dependent thermal conductivity 

In the following we compute the quasiparticle density of states and the optical conductivity 
for the f-wave gap and for the multigap model in the Born limit at T = OK, assuming, that 
both inter- and intraband scattering of impurities are important, which was predicted to have 
strong effect on the superconducting state [22] ■ The density of states differs significantly in 
these two models, producing distinct peaks from the active and passive bands in the ZR model, 
while only one single peak is found in the f-wave case. In the optical conductivity, beyond the 
intraband resonance peaks another extra feature is found coming from interband scattering. 
On the other hand, the single broad bump in the f-wave case can clearly be distinguished from 
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the three peak structure of ZR's model, hence both the single particle density of states and 
the optical conductivity manifest evident signatures of the corresponding model. 

Density of states. - As already pointed out by Agterberg [22], the effect of intcrband 
impurity scattering plays an important role in the multigap model, similarly to the interband 
scattering of Cooper pairs, which induced A2(k) in the passive bands. However, in the f-wave 
case, we deal with a simple one band model, therefore only intraband impurity scattering is 
considered. 

The effect of weak impurities (Born limit) is incorporated in the quasiparticle Green's 
function of the multigap model by renormalizing the frequency: 



v 1 ^i V v 1 - u l J v l - us 



A 2 m= w + r u - it —^ k\ -^=— ) +r 12 ^==, (4) 



where u\ and m 2 are the dimensionless frequencies renormalized by A 2 and A p in the passive 
(A 2 ) and active (A p ) channels, respectively, Tji = niirNi\Uji\ 2 = 1,2), where the total 
density of states is split as Ni : 7V 2 = 0.43 : 0.57 between the passive a + (3 and the active 7 
bands from de Haas- van Alphen measurement |23| . V\\ and U22 are the intraband, U\2 = U21 
are the interband scattering matrix elements, n% is the impurity concentration, K{z) is the 
complete elliptic integral of the first kind. As opposed to this, for the f-wave model, the 
frequency is renormalized as 

^ « ,k(*), (6) 



Af A/7T VI -u 2 Wl -u 2 

where T = ni7riVo|£/| 2 , ^0 = ^1 + N2 is the total density of states. At the Fermi energy 
(oj — 0), the frequencies are obtained as 

Ul = 4iexp (f A2 r ll( A p X) + r 12 r 21 ) ' < 7a > 

U2 "a~T7i — r VTr r Ul ' (7b) 
1 ul 1 — J- 22J + J- 121 21 

u = 4iexp(^— (7c) 



, 2r , 

for small T's. The quasiparticle density of states (DOS) is obtained as 

AWigapM = Im ( Nl ^p==K l^=L=) + N2^=) , (8) 

N f (w) = N Q lm^L=K (^=L=) (9) 

for the multigap and f-wave model, respectively. In the Born limit, the importance of interband 
scattering readily follows from eq. Q. Without interband processes, the two bands decouple, 
and the p-wave gap behaves like an s-wave gap in the presence of nonmagnetic impurities: 
it remains gapped for small impurity concentrations, and reaches the gapless region only for 
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high concentrations or strong scatterers. In the density of states, peaks are expected at the 
gap maxima of the two bands, namely at to = A 2 and A p . As soon as interband scattering 
is present, the p-wave channel becomes gapless for arbitrary concentration and strength, as 
is seen from eq. In figs. |5]and|21 the density of states is plotted for U\\ = U 22 = U, 
T = O.OOlAp and for weak (T 12 N /N 2 = 0.000lA p ) and strong (T 12 N /N 2 = 0.1A P ) interband 
scattering. Here we used A2 = 0.3A p following refs. |19H21| . which was found to describe 
very well the experimentally obtained specific heat and penetration depth. Also for the f- 
wave model, the weak-coupling relation Af = 1.2lA p was assumed. In the weak interband 
scattering limit, the DOS in the 7 band increases very slowly, while for strong impurities, even 
a small peak develops at u> = A 2 in the p-wave DOS, due to the strong interplay of the two 
bands. For the f-wave model, only one peak is present at lo = Aj, henceforth clear indications 
of either the multigap or the f-wave model should be observed experimentally. 




0.2 0.4 0.6 0.8 1 1.2 14 1.6 1.8 2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



U! I Ap OJ I Ap 

Fig. 2 Fig. 3 

Fig. 2 - The density of states is shown for weak interband scattering for the multigap model: in 
the passive band (dashed-dotted line), active band (dashed line) and the full DOS (solid line). The 
corresponding f-wave DOS (dotted line) exhibits one single peak at to — 1.21A P . 

Fig. 3 - The density of states is plotted for strong interband scattering. The peak in the p-wave 
channel (dashed line) at lo = 0.3A P indicates the influence of these kind of processes. 

In fig. ^ the DOS is shown for medium strength inter- and intraband scatterers: Uji = U , 
r = O.lAp. The peaks at the gap maxima are broadened into big bumps and increasing 
density of states is observable in the p-wave channel below the gap maximum. In fig. the 
DOS of both models are compared, as intraband scattering increases. 

Optical conductivity. - In the following we consider the optical conductivity of the multi- 
gap and f-wave model at T = OK. In the former case, as a lowest order approximation, each 
band contributes separately to the optical response, which is obtained as 

Rea(u) = CiRecri(cj) + ( 2 Rca 2 (u;), (10) 

where £i = (ni/mi)/(n/m) = 0.64, £ 2 = (n 2 / m 2 ) / (n / m) — 0.36, which measures the con- 
tribution of each band to the conductivity • At T = OK, the conductivities are obtained 
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Fig. 4 
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Fig. 4 - The density of states is shown for medium strength impurities for the multigap model: in 
the passive band (dashed-dotted line), active band (dashed line) and the full DOS (solid line). The 
corresponding f-wave DOS (dotted line) exhibits one single peak at u — 1.21A P . 

Fig. 5 - The density of states is shown for the multigap (solid line) and f-wave model (dotted line) 
for Uu = U 22 = U,r = O.OOlAp, 0.01A P and 0.1A P , r 12 N /N 2 is fixed to 0.01A P . 



as 



Re<77(w) 



2uim 



Re 



^Fi[ui(lo - x), -ui(x)] - Fi[ui{uo - x), ~ui(x)]^ dx, 



(11) 



where [2HI2H1 
Fr{u,u') = 



7rA 2 a'-u 1 yJY 
F 2 (u,u') = 



,'2 



-,K 



,'2 



-.K 



1 



1 



A p u' 



For the f-wave model, the conductivity reads as 



vT 



Recr(w) 



2ujm 



Re 



— x), — u(x)] — F\[u{uj — x), — 



dx, 



(12) 
(13) 

(14) 



and A 2 in Fi(u, u') has to be replaced by Af. 

^From the density of states of the multigap model, one would expect peaks in the conduc- 
tivity at u> — 2A2, 2A p , and if interband scattering is strong enough, at u> = A p + A 2 . These 
features can bee seen in figs. and \7\ Here the intraband scattering parameter was fixed 
as Uu = U22, FiiNq/Ni = O.OOf Ap and O.Of A p , respectively, and the interband scattering 
parameter varied from weak through medium to strong as r 21 = O.lFn, Tn and lOFn. Along 
this line, as r 21 enhances, the peak corresponding to interband excitations at u) — A p + A 2 
becomes dominant, overwhelming the intraband processes. 
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Fig. 7 - The optical conductivity of the multigap model is plotted for Uu = U22, TuNo/Ni = 0.01A P , 
the interband scattering parameter changes as F2i = O.irn, Tn and lOFn, with increasing peak at 
u> — 1.3A P . 




Fig. 8 - The optical conductivity of the f-wave model is shown for TuNo/Ni = 0.001A P (solid line, 
enlarged by 50), 0.01A P (dashed line, enlarged by 12) and 0.1 A p (dashed-dotted line). 

Fig. 9 - The optical conductivity of the multigap model is shown for r^iVo/A^ = 0.01A P , FuNa/Ni = 
O.OOlAp (solid line, enlarged by 50), 0.01A P (dashed line, enlarged by 12) and 0.1 A p (dashed-dotted 
line) . 



In the f-wave case, the conductivity is identical to those of a d-wave superconductor, and 
exhibits a small peak at w — as the impurities get stronger, and at the same time a rather 
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broad bump is developed at w = 1.5 — 2A p as is seen in fig. [S] It is remarkable also, that 
in spite of the different scattering rates changing within two orders of magnitude, the curves 
look rather similar to each other, and can easily be distinguished from the multigap model. 
For the multigap model, as intraband scattering becomes dominant, the conductivity of the 
passive band rules over the other transitions, since the same impurities decrease the superfluid 
density much faster in this band due to the smaller gap (see fig. 

Conclusion. We have calculated the density of states and the optical conductivity 
for ZR's multigap model and for the f-wave case in the Born scattering limit. The effect of 
interband scattering produces new structures in the optical response of the multigap model: 
resonant peaks can be observed at u> = 2A 2 , 2A p , and due to interband transitions at A2 + A p . 
Therefore when the experimental difficulty to perform the optical measurement below T = 
O.IK is overcome, this will provide a unique way to explore the validity of models for the 
superconductivity in Si'2Ru04. In a forthcoming paper we shall discuss the effect of resonant 
scatterers, where we also expect clear differences between the multigap and f-wave model. 
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